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^ Abstract 

We use lattice simulations to study the single-site version of SU (N) lattice gauge theory with 
two flavors of Wilson-Dirac fermions in the adjoint representation, a theory whose large volume 
correspondent is expected to be conformal or nearly conformal. Working with as large as 53, we 
map out the phase diagram in the plane of bare 't Hooft coupling, g'^N, and of the lattice quark mass, 
" am, and look for the region where the center symmetry of the theory is intact. In this region one 

expects the large- equivalence of the single site and infinite volume theories to be valid. As for the 
Nf = 1 case (see Phys. Rev. D 80: 065031), we find that the center-symmetric region is large and 
includes both light fermion masses and masses at the cutoff scale. We study the A^-dependence of the 
width of this region and find strong evidence that it remains of finite width as A^ — )• oo. Simulating 
with couplings as small as g'^N = 0.005, we find that the width shrinks slowly with decreasing g'^N, 
at a rate consistent with analytic arguments. Within the center-symmetric region our results for 
the phase structure, when extrapolated to N = oo, apply also for the large volume theory, which is 
minimal walking technicolor at A^ = oo. We find a first-order transition as a function of am for all 
values of b, which we argue favors that the theory is confining in the infrared. Finally, we measure 
the eigenvalue densities of the Wilson-Dirac operator and its hermitian version, and use large Wilson 
loops to study the utility of reduction for extracting physical observables. 
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I. INTRODUCTION 



There has been a recent revival of interest in the possibihty of using complete volume 
reduction for the infinite N (number of colors) limit of QCD and QCD-like theories. If this 
reduction holds, then the theory, defined nonperturbatively on a lattice, gives predictions 
that, at infinite A^, are independent of the number of sites. Specifically this means that the 
theory defined on a single site, or a small, fixed number of sites, is large- equivalent to the 
corresponding infinite volume theory with the same bare parameters 

Reduction to a single-site has long been known to fail for the pure gauge theory (and thus 
also for QCD in 't Hooft's large- iV limit, since quark contributions are suppressed by in 
this limit) l2-4|. This failure is due to the breakdown of one of the conditions needed for a 

n 

largCj^A^ orbifold equivalence between the single-site and large volume theories (see Refs. p] 
and [6]). This key condition is that the Zjf center symmetry of the single-site theory must 
be unbroken. This symmetry breaks spontaneousljll] in the single-site pure gauge theory (the 



Eguchi-Kawai [EK] model [l!]). This is expected from perturbation theory (PT), where, at 
leading order, the effective potential for eigenvalues of the holonomy around the compact 
direction (the Polyakov loop) leads to attraction and thus clumping j2|, |3| . Several years ago, 
it was realized that the addition of massless fermions that reside in the adjoint representation 
and that have periodic boundary conditions in the compact directions leads (in perturbation 
theory) to a repulsion between eigenvalues, which in turn leads to a uniform distribution of 
these eigenvalues In this case the center symmetry is unbroken and reduction holds. 

Two of us have previously investigated the single-site theory with a single Dirac adjoint 
fermion (discretized using Wilson fermions), finding that, for small to rather large values of 
the inverse 't Hooft coupling, b = l/{g^N) E [0, 1], there is a large range of values of the quark 
mass for which the center symmetry appears to be unbroken j^. This result was unanticipated 
because several leading-order perturbative calculations, done with a single compact direction 
(or a single site along only one Euclidean direction) show that the symmetry breaks if the 
physical mass m exceeds a value of 0(1/ aN) Ulnl. As m increases from ~ 1/aN to ~ 1/a, 



^ Strictly speaking, the symmetry is spontaneously broken only for N ^ oo. In practice, however, effective 

spontaneous symmetry breaking is seen in simulations at finite but large values of N, and we use the 

terminology of phase transitions throughout this article. 
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the perturbative calculations indicate that the eigenvalue density of the link in the short 
direction will form a number of clumps, starting with 0{N) clumps at very small masses, and 



decreasing to a single clump at infinite mass. The results of Ref. 



(which used up to 



15) do not follow this pattern, instead finding no clumping for masses up to of 0{l/a) for all 



values o 
Refs. UJ 



A^. A semi-quantitative understanding of these results has recently been given in 



12] • With more than two compact dimensions, the fluctuations in the eigenvalues can 



overwhelm the tendency to clump, and this happens up to masses of 0(1/ a). The numerical 



results for the phase diagram obtained in Ref. 



Ref, 



have also been checked, and extended, in 



We also note that simulations with Nf = 1/2, 1 and 2 massless overlap adioint 

nn 

fermions also find no center-symmetry breaking, at least at large b [131. Il4j|. 

In the present paper we extend our investigations to Nf = 2. The main motivation for 
doing so is that the corresponding infinite volume theory is expected to be nearly conformal, 
and thus a candidate "walking technicolor" modelj^ Indeed the theory with two colors is the 
theory with the smallest field content that lies close to the conformal window, and has been 
dubbed the "minimal walking technicolor" (MWT) model. If reduction holds, then we should 
be able to study a close relative of MWT, i.e. the theory with A^ = oo. One naively expects 
only a weak dependence on A^, because the number of both gluonic and fermionic degrees of 
freedom scale as A^^. 

We also note that the Nf = 2 AEK (Adjoint Eguchi-Kawai) model is related by a combi- 
nation of orbifold and orientifold equivalences to the QCD-like theory with 2Nf = 4 fermions 



in the two index anti-symmetric (AS) irrep 



This theory in turn is the large- A^ limit 



of QCD with 2Nf quarks, but with the limit taken with the quarks in the AS irrep (which is 
equivalent to the anti-fundamental for A^ = 3). This is the Corrigan-Ramond large- A^ limit of 
QCD 18], which differs from the 't Hooft limit in having fermion loopsQ 

Our main effort herein is to determine the phase structure of the Nf = 2 AEK model. 
To do so we have upgraded our simulation algorithm from a Metropolis algorithm, with CPU 
scaling as A^^~^, to a Hybrid Monte-Carlo (HMC) algorithm, for which we find CPU oc A^^~^'^. 



15j and 



16|. 



^ For recent reviews of technicolor models on the lattice, see Refs. 
^ This equivalence holds only in charge-conjugation even subsectors. 

^ Note that here taking the Corrigan-Ramond limit moves one from a theory which is not close to the conformal 
window (4 quarks in the fundamental irrep) to one that is (4 quarks in the AS irrep) , suggesting that 1 /N 



corrections are probably large, at least in this respegt. 



This allows us to reach much larger values of A^, and to improve the statistics. Together, these 
advances allow us to study the nature of the symmetry breaking in more detail than in the 



Nf = 1 study 
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12|. 



allowing us to compare with the theoretical expectations of Refs. 
Our main result is that we find the phase diagram to be qualitatively similar to that for 
Nf = 1, with center symmetry remaining unbroken for masses up to 0(l/a). Specifically, our 
evidence suggests that, at fixed coupling, although the range of masses for which the symmetry 
is unbroken shrinks somewhat as increases, it remains of width ~ 1/a as — )■ oo. Our 
strongest evidence for this is at 6 = 1, but our results suggest that this holds for 6 = — 200, 
i.e. for the entire range of coupling that one could possibly be interested in. Thus our results 
suggest that one can use adjoint fermions of almost any mass to "stabilize" reduction. This 
is only expected to fail in the extreme weak coupling limit {b — ?■ oo). 

This is an encouraging result, and so we have made the first steps in trying to see if 
reduction can be used to obtain results for physical quantities. The key question is how large 
a value of N is needed so that the physical contributions are larger than those from 1/A^ 
effects. We have investigated this by studying the large extrapolation of the plaquette, by 
calculating the spectrum of the Wilson-Dirac operator and its hermitian counterpart, and by 
calculating large Wilson loops in order to see if we can extract the heavy-quark potential. 



Work along similar lines has recently been reported in Ref. 19|. These authors simulate 
the theory with two adjoint Wilson fermions on a 2"^ lattice, and use A^ = 2 — 6. They report 
evidence that, for 6 = 2, there is a region of quark masses around the putative critical value, 
including quarks of masses ~ 1/a, for which the center symmetry is unbroken. 

This paper is organized as follows. In the following section we describe the AEK model 
and the properties of the large- volume theory to which it would be equivalent were reduction 
to hold. In Sec. IIIII we describe the algorithm that we use, and show some results concerning 
its performance. Section HVl is the core of the paper, in which we use our numerical results to 
determine the phase diagram of the AEK model. We then, in Sec. |Vl present first results for 
"observables" — the spectra of the Wilson-Dirac operator and its hermitian counterpart, and 
large Wilson loops. We close in Sec. I VI I with a summary and a discussion of the outlook for 
future work. An appendix describes models for eigenvalues of the Wilson-Dirac operator that 
are used in Sec. |Vl 
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II. THE AEK MODEL AND ITS PUTATIVE LARGE- VOLUME EQUIVALENT 



The partition function of the single-site theory is 

Zaek = / Y\DU^Di/jDiI) exp S'gaugc + D^'4^j , (2.1) 

V J 

where the four C/^ are SU(N) matrices, while ■0j and ipj are Grassmann Dirac variables of 
flavor j, living in the adjoint representation of SU{N). We use the Wilson gauge action 

W ^2NbJ2 ReTr U^l^ , (2.2) 

where UJ^]^'^ is the product of links around the plaquette in the ii, v plane, and h is the inverse 
't Hooft coupling. 

We also use Wilson's lattice Dirac operator 



Di^ — \ — K 
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(2.4) 



where f/^"^' is the adjoint representative of U^^ and k is the usual hopping parameter, related 
to the bare quark mass by 

rno = ^ - 4 . (2.5) 

Periodic boundary conditions on both gauge and fermion fields have been built into the form 

of -DvK- Throughout this paper we set A^j = 2. 

The theory has a center symmetry, under which 

V^^z''-U^, (2.6) 

where z = exp(2z7r/A^) and Q < < N are integers. Note that is invariant under 
this transformation, so that the fermion action is also invariant. There is also the single-site 
version of the gauge symmetry 

nu^n^ [n e su{n)] . (2.7) 



Finally, there is an SO (4) flavor symmetry, most easily seen by writing the action in terms of 
four Majorana fieldslfl 

If reduction holds, this single-site theory is equivalent, when — )■ oo, to a theory that 
has any number of lattice sites in each of the periodic directions /i = 1,2,3,4, including 
the case of = oo. The action of the > 1 theories has the same form except that U, 
ip, and ip, are now fields having a site index, S'gauge contains a sum over the position of the 
plaquettes, and connects fermion fields at adjacent sites. We stress that an important 
feature of reduction is that it relates the single-site and infinite-volume theories having the 
same bare parameters, b and k. 

We recall some important properties of the infinite-volume theory, since these will be 
inherited by the single-site theory if reduction holds. First, the theory is asymptotically free — 
Nf > 11/4 fermions are required to change the sign of the first coefficient of the /3-function. 
Second, although the bare quark mass vanishes when k, = 1/8, this critical value of k, is 
additively renormalized because Wilson fermions do not preserve chiral symmetry. The critical 
value is shifted to Kc{b) > 1/8, and the physical quark mass becomes 

Here we have introduced the lattice spacing, a, which can be determined, in principle, by fixing 
the value for a physical scale, such as a particle mass. Since the theory is asymptotically free 
at short distances, one approaches the continuum limit (a — ?■ 0) by sending 6 — oo, and in 
this limit Kc{b) -> 1/8. 

The nature of this critical line depends on the infrared behavior of the theory. One pos- 
sibility is that the theory lies below the conformal window, so that chiral symmetry is spon- 
taneously broken, much as in QCD. Then, for n near one can study the long-distance 
behavior and vacuum structure of the lattice theory using chiral perturbation theory (ChPT). 
In particular, close to the continuum limit, one can use a modified ChPT which includes dis- 
cretization effects 20|]. For adjoint fermions, the symmetry breaking pattern differs from that 
in QCD, and is S'[/(4) — )■ SO {4). The required generalization of the analysis of Ref. 20[ has 



^ In the continuum this symmetry becomes an S'C/(4) symmetry, but, with our choice of fermion discretization, 
only its 5*0(4) subgroup remains an exact symmetry of Eq. (|2.ip ). 
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been given in Ref. 21j. One finds that, as in QCD, there are two possible scenarios: either 
there is a first-order transition hne, at which the degenerate pseudo-Goldstone "pions" attain 
their minimal, non-zero mass, or there are two second-order lines, along which the pions are 
massless, and between which there is an Aoki-phase 22|. Within the Aoki-phase, the SO {4) 
vector symmetry is broken]^ The width of the Aoki-phase is oc a^, and thus shrinks rapidly 
as one approaches the continuum limit. 

A different possibility for the critical line arises if the massless theory is conformal in the 
infrared, i.e. if there is an infrared fixed point. There is growing numerical evidence that this 
is the situation in the N = 2 theory. For this theory, the simulations of Refs. {25 1 and 26 1 
map out parts of the phase diagram. In particular, there is single, second-order transition line 
emanating from {b = 00, k = 1/8), while for 6 ^ 1/4 the line becomes a first-order transition. 
(A similar picture holds for an improved fermion action, but in this case the second-order 
line extends to stronger coupling 27|].) This is not established definitively, and also does not 
directly apply to the N = 00 theory that we are interested in. Nevertheless, this possibility 
provides a quite different phase diagram than that which applies when one is outside the 
conformal window. One of our aims is to see which possibility holds at = 00 (assuming 
that reduction holds). 



III. SIMULATION ALGORITHM 



We simulate the single-site theory using the hybrid Monte-Carlo (HMC) algorithm 
Integrating out the fermions leads to 

Zaek = / n ^f^/^e^'^"- det{Dwf • (3.1) 

As usual, 75 hermiticity implies that det{D\y) is real, so we can write 

det{Dwf = det{DwDl,) = det{Q^) , (3.2) 

where Q = Dw^5 = is the Hermitian Wilson-Dirac operator. Since has positive 
eigenvalues, we can represent its determinant using pseudofermions. Introducing momenta 



^ We note for completeness that Refs. 
"first-order scenario." 



24| have recently raised concerns about the consistency of the 
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conjugate to the link variables, we end up with the HMC Hamiltonian 



H=IY.^t {PD - 2Nb ReTrt/P^ + <P^Q-'<P . 



(3.3) 




The are traceless hermitian N x N matrices, while the pseudofermion is complex, lives in 
the adjoint representation of SU{N), and has an implicit Dirac index. It thus has 4(A^^ — 1) 
complex components. 

In practice, we represent in color space as a traceless bifundamental, traceless 



In this way we do not need to explicitly construct f/^'^-'l^ In fact, one could, in principle, keep 
the trace of 0, since the singlet field that it represents has no impact on the dynamics. In 
particular, one can show that, in exact arithmetic, the molecular dynamics (MD) trajectories 
that are followed are identical with or without tr included, as is the change in H. We find, 
however, that the number of CG iterations required for a given accuracy is larger if tr0 is 
included, presumably because one has to do some work to find the solution for the singlet 
part. Thus we always set tr0 = 0. 

Our implementation of the HMC algorithm is standard. We invert Q"^ using the conjugate 
gradient (CG) algorithm, with a weaker stopping criterion during the MD evolution than for 
the accept-reject step. We require that the residue, r = b — Q^x, with b the source, satisfies 
|rp/|6^| < 10^^ during MD evolution, finding that any further increase of the cut-off leads to 
a drop in the acceptance. For the accept-reject step we use < 10~^^, which makes the 

error in AH negligible. Our CG always starts from a vanishing guess, Xq = 0, which assures 
reversibility of the trajectory. We use trajectories of unit length, and adjust the step size to 
attain acceptances of 0.6 — 0.85. 

In deriving the gluonic force, one must account for the fact that each link appears twice in 
each plaquette. Nevertheless, the final result has the standard large-volume form: 



N X N matrix, on which U^'^^ acts as 




(3.4) 




(3.5) 



^ We thank Simon Catterall for stressing this point to us. 
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The fermionic force is 

-(75 - 1,1,U [U.i'pUlxi - xlU.i^pUl] I + h.c. (3.6) 

where a and /3 are Dirac indices, x = Q '^'Pi and if) = Qx- Both forces maintain the traceless- 
ness of Pfj,. 

We now discuss the scahng of CPU time with A^, which is a key factor in determining how 
large one can take N. The core operation — multiphcation of N x N matrices — scales as N^. 
The use of the bifundamental form of U^, rather than the adjoint, is crucial here, reducing the 



scaling from A^^ to A^^ , as pointed out in Ref . 



191 ] ■ The next contribution to the overall scaling 



comes from the number of CG iterations, Ncg- This turns out to depend on the proximity 
to the critical line. An example is shown in Fig. [H for the stopping criteria given above. 
Away from the critical line, Ncg is independent of A^, while near the line it grows roughly like 
A^^/^. The third ingredient is the inverse step size, or equivalently the number of MD steps 
(Nmd) per trajectory (for a given acceptance rate). We find that, to good approximation, this 
grows linearly with A^. Thus, for trajectories of unit length, CPU time scales as ~ A^^ away 
from the critical line, and roughly as ~ A^^'^ near to the line. Both scalings are considerable 
improvements over that for the Metropolis algorithm used in Ref. Q|, which is A^^ for each 
SU(2) subgroup update and A^*^ for an entire update. On the other hand, our scaling is not 



as good as the estimate of A^^'^ given in Ref. 19|, which assumed Nmd oc N^^"^ and that Ncg 



is independent of A^, and explicitly excluded the possible effects of critical slowing down. 

We have done both horizontal (fixed h) and vertical (fixed k) scans in the b — k plane, 
studying the ranges k, = — 0.6 and b = 0.05 — 200, although our main focus has been on the 
smaller ranges k = — 0.26 and b = 0.35 — 1.0. We use A^ = 10 — 30 in these scans. Rather 
than quote a complete set of run parameters we give a representative example. We have, for 
b = 0.35, 0.5, 0.75 and 1.0, used 27 values of k (0.0 — 0.26 in steps of 0.1). At each n, we start 
from the "configuration" output from the previous value, thermalize for 500 trajectories, and 
then run for 7500 (A^ = 10), 5000 (A^ = 16), 2000 (A^ = 23) or 1000 (A^ = 30) trajectories 
during which we make measurements every 5 trajectories and store the configurations every 
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FIG. 1. Average number of CG iterations in the MD updates for various as a function of k at 
b = 1.0. Results are from UP scans. 

50. Each scan is done in both directions — the UP and DOWN scans denoting increasing or 
decreasing parameter values (either n or b). To give an example of the CPU time required, 
the b = 1.0 UP scan took 33, 155, 342 and 618 CPU-hours of a single core on 3.0 GHz Intel 
Xeon processor, for = 10, 16, 23 and 30, respectively. Our simulations have been done on 
local workstations and using up to 32 CPU cores on a computing cluster. 

We have also done longer runs at several points in the b — k plane, in which we have gone 
up to A^ = 53. Details of these runs will be given below. 



IV. PHASE DIAGRAM OF THE Nf = 2 AEK MODEL 



In this section we present our main results, from which we deduce the phase diagram 
sketched in Fig. |2J The most important conclusion is that there is a "funnel" in which the 
center symmetry is unbroken, on either side of the first-order transition which we identify with 
Kc- The diagram is qualitatively similar to that found for Nf = 1 8|. 
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0.04 0.08 0.12 0.16 0.20 0.24 

FIG. 2. Sketch of phase diagram for the A'^^ = 2 AEK model in the k — b plane for N ~ 30. Note 
that the k = axis is the EK model. The positions of phase boundaries are approximate, and 
depend somewhat on A^. The shaded region at Kc indicates the uncertainty in the position of what 
appears to be a first-order transition due to hysteresis. Within each region we note the subset of 
the Zff center symmetry that is unbroken, with Zi indicating complete breakdown. The center 
symmetry is unbroken in the hysteresis region. The detailed symmetry-breaking pattern for large k 
is representative, and depends to some extent on N. For further discussion, see text. 

A. Measured quantities 

To study the gross features of the phase diagram, we calculate the average of the plaquette, 
Up, defined by 

^r^^J2Tr{K'r)- (4.1) 
10 



To study center symmetry breaking, we consider general "open loops" : 

= 1-tT U^^ U^' U^' U2\ with n^ = 0,±l,±2,... (4.2) 

where = ?7^". These loops transform non-trivially under the center symmetry, unless 
all four TL, are integer multiples of N. They are thus order parameters for center-symmetry 
breakingjj The simplest choices, on which we focus, are the four Polyakov loops, = -^tr t/^ 
and the 12 corner variables, M^i, = -^trf/^f/^ and M^^^^ = j^tiUfJJl, with ii^v. 

As in the quenched Eguchi-Kawai model, the corner variables turn out to be particularly 



29). 



We illustrate this with 



useful because they are sensitive to partial symmetry breaking 
simple examples. First, if f/^ = 1 for all /x, then all the Polyakov loops and corner variables 
are unity. This corresponds to complete breaking of the center symmetry. If instead 

V/i: [/^ = diag 1,...,1, -1,...,-1 (4.3) 

\iV/2 entries N jl entries / 

(where we have assumed that is divisible by 4, so that det f/^ = 1), then the Polyakov loops 
vanish, while all corner variables are unity. In this case, there is a unbroken subgroup: both 
(P^) and (M^ -i-j,) are invariant under the Zi subgroup of generated by 

V/i : U^^ ±U^. (4.4) 

The themselves are invariant under the combination of Eq. fl4.4p and a gauge transforma- 
tion, the latter being the similarity transformation which interchanges the first A^/2 diagonal 
entries with the second A^/2 entries. 

Such partial symmetry breaking can be discussed in a gauge invariant way by considering 
the eigenvalues of link matrices. For each link we can write 

U, = W^K^Wl , e SU{N) , (4.5) 

with containing the eigenvalues: 

A, = dmg(e<,e''i,...e<) . (4.6) 



We have used —5 < < 5 to keep the quantity of data manageable. This has the disadvantage that the 
traces are then insensitive to symmetry breaking such as Zfj — > Ziq. Histograms of hnk eigenvalues, to be 
discussed below, are, however, sensitive to such symmetry breaking. 
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Gauge transformations can permute the eigenvalues, but not change their values. Center- 
symmetry transformations change the eigenvalues by a uniform translation: 0^ — > + 
27Tn^/N. Thus a direct way of looking for certain symmetry breaking schemes, and under- 
standing their nature, is to look at the distributions of the 9^. For example, unbroken center 
symmetry implies a distribution which is invariant under translations by 27m/ N. Partial sym- 
metry breaking occurs when a subgroup of such translations is unbroken. In the first example 
above, the eigenvalues are all clumped, and all translation symmetries are broken. In the 
second example, the eigenvalues form two clumps, and translation by n remains a symmetry. 

We use the link eigenvalues in Sec. lIVEt plotting histograms and considering the correla- 
tions between links in different directions. 



B. Scans at moderate coupling (6 < 1) 

In this section we use scans of the plaquette, Polyakov loops and corner variables to map 
out the gross features of the phase diagram. 

The most interesting values of b are roughly 0.35 — 1.0; this was the range studied in the 



Nf = 1 model [s]. For = 3, this corresponds to /3 = Q/g"^ = 6.3 — 18, a range running 
from couplings similar to those used in large-volume simulations to very weak coupling. We 
have made detailed scans at 6 = 0.35, 0.5, 0.75 and 1.0. That at 6 = 0.35 shows a great 
deal of structure that is hard to analyze (including large hysteresis and the influence of a 
bulk transition), while that at 6 = 0.75 interpolates between the results at 6 = 0.5 and 1.0. 
Thus we show, in Figs. |3] and H] respectively, scans of the plaquette at 6 = 1.0 and 0.5. We 
have simulated with N = 10, 16, 23 and 30, but, for the sake of clarity, show results only for 
= 16 and 30. We also show, in the central region, an approximate estimate of the result 
at = oo, obtained by fitting results at the four values of A^ (or more values, if available) to 
Co + Ci/N + C2/N^- Such fits will be discussed in Sec. IIV CI 

The results at 6 = 1 show three main features: (i) a change in slope at k, ^ 0.02 (and 
possibly another at k ^ 0.05), (ii) a jump at k = 0.13 — 0.14, and (iii) a transition region 
at K ~ 0.2. These correspond on the phase diagram of Fig. [2] to (i) the transition region 
from center-symmetry broken phases to the unbroken central region, (ii) the transition line at 
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FIG. 3. Scans (both UP and DOWN) of the average plaquette at b = 1.0 for = 16 and 30. The 
results of an extrapolation to = oo (described in the text) are shown in the central region. 
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FIG. 4. As for Fig. H but for b = 0.5. 
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FIG. 5. Scans of the absolute values of corner variables for b = 1.0, for = 16 and 30. Results for 
the 12 independent |M^jy|'s are shown separately. 

~ 0.13, and (ill) the transition to the region of multiple broken phases for large k. We focus 
first on the central feature, presenting our evidence concerning symmetry-breaking below. An 
important issue is whether the jump at finite survives as a first-order transition at = oo. 
Our extrapolations suggest that it does: although the jump in the plaquette decreases with 
A^, it appears to remain finite at N = oo. 

The conclusion of a first-order transition is clearer in the results at 6 = 0.5. These show 
the same qualitative features as for 6=1, but the jump in the plaquette is larger, and there 
is hysteresis for A^ = 16 and 3o|^ 

To study center symmetry breaking, we first use scans of the absolute values of Polyakov 
loops and the corner variables. Both should vanish as A^ — )■ oo if the symmetry is unbroken. 
The corner variables are more informative and we show an example, for b = 1.0, in Fig. |3 
^ We note in passing a peculiar phenomenon we have seen for smaller values of N in the range we consider. 
At N = 10, the UP scans for b = 0.5 show, in the hysteresis region, points having average plaquettes 
with non-vanishing imaginary parts. This breaks the charge conjugation symmetry of the theory, and is 
reminiscent of results found in the twisted EK model 30| . It is because of these points that we do not have 



an extrapolated result for the UP scans at k = 0.15 — 0.17 in Fig. We suspect that this occurs only in 
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metastable phases, and view it as a sign of the complicated vacuum structure of the single-site theory. 



Results for 0.35 ^ b < 1.0 are qualitatively similar. 

For both small and large k, k ^ 0.05 and k ^ 0.19, the corner variables indicate that the 
center symmetry is broken. The nature of this breaking is clarified by the Polyakov loops, 
\Pfj,\, whose plots we do not show for the sake of brevity. For n ^ 0.02 we find |P^| to be 
non- vanishing as — )■ oo, indicating that the center symmetry is completely broken. This is 
the Zi phase shown in Fig. [2l For 0.02 ^ n ^ 0.05 and k = 0.19 — 0.20, however, |P^| are 
consistent with zero at = oo, indicating only partial symmetry breaking. The nature of 
this partial breaking can be elucidated using the distributions of and M^,y in the complex 
plane, and using histograms of link eigenvalues. Some examples of the latter will be shown in 
Sec. HVEl 

A key issue for reduction is the realization of center symmetry in the central funnel, 0.05 ^ 
K, ^ 0.18. For the values of A^ used in the scans, we find no indication of symmetry breaking. 
Our evidence is as follows. Scatter plots of {M^y) and (P^) in the complex plane show a single 
distribution centered around the origin, with averages consistent with zero. Similarly, the 
higher-order traces, Kn-, which we have evaluated at several positions inside the funnel, are 
all consistent with zero. Finally, histograms of link eigenvalues, examples of which are shown 
in Sec. IIVEI are also consistent with the absence of symmetry breaking. 

The other key question is whether the funnel remains of finite width as A^ — ?■ oo. We can 
see from Fig. |5] that the funnel narrows with increasing A^. In particular, the lower edge of 
the funnel, which we call kj, increases from hif ^ 0.03 at A^ = 10 to ~ 0.05 at A^ = 30. 
To study this question further requires larger values of A^, and we defer consideration until 
Sec. UVDl 

To complete the study of the phase diagram we have done several vertical scans with 
A^ = 16, 23 and 30. We show an example of the results in Fig. [6l which displays the (|M^j,|) 
for A^ = 23 and k < For all i^i there is no indication of symmetry breaking at strong 
coupling, h ^ 0.3, just as for the EK model. There are possible transitions, however, as we 
increase h. For example, at k = 0.02 we see two transitions: one at 6 ~ 0.3 and a second at 
h ~ 0.55. The first is from a center symmetric phase to one in which both Polyakov loops 
and corner variables are non-vanishing, consistent with complete symmetry breakdown. The 
second is to a phase with large corner variables and smaller Polyakov loops, which we interpret 
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FIG. 6. Vertical scans of the absolute values of corner variables with = 23 and for k = 0.02, 0.04, 
0.06 and 0.12. For k = 0.02 both UP and DOWN runs are shown. 

as a partially broken phase. This is the same "Z2" phase apparent in Fig.[5]for k, = 0.02 — 0.04. 

For all other k there is no hysteresis, so we show only UP scans. At k, = 0.04, we see 
only a single transition, at 6 ^ 0.3, and this is directly to a partially broken phase with only 
{\Mfj_i,\) non-zero. For higher n, however, the symmetry is unbroken on both sides of the jump 
at 6 ^ 0.3, and we interpret this as a bulk transition. It is unclear, however, whether this 
corresponds to a phase transition or a crossover as — i- 00. 

Vertical runs at higher values of k fill in gaps left by the horizontal scans, and are part of 
the input which leads to the phase diagram of Fig. [2l For the sake of brevity, however, we do 
not show any plots here. 



C. Scaling of the plaquette, (|-P^P), {{M, 



|2\ 



In order to study the key question of whether the symmetry unbroken funnel remains as 
N ^ 00, we have extended the calculations to larger values of N at several several values of 
b and k. 
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We begin by looking at the average plaquette. If reduction holds, then, away from Kc, 
the single-site theory is equivalent at large to a large-volume lattice theory with quarks 
having physical masses of 0(1/ a). The long-distance physics of such a theory is that of a 
pure-gauge theory with action modified from the pure Wilson form by fermionic effects. If 
K is much smaller than Kc, these modifications should be small, since they are proportional 
to powers of the small quantity k (as can be seen using the hopping parameter expansion). 
The large-volume theory is thus close to the pure-gauge theory with Wilson action. We can 
therefore make the semi-quantitative prediction that, near the lower boundary of the funnel, 
Kf, the average plaquette should lie close to the large- volume, pure-gauge (Wilson action) 
value, but depart from that value as one approaches k^. On the other side of the transition, 
however, we do not expect the plaquette to tend to this same value as k — Kc increases. This 
is because k is now larger, so the action differs more significantly from the pure-gauge Wilson 
form. 

Figures [3] and m show that the plaquette has considerable dependence on A^, with the slope 
of this dependence varying with n. We show in Fig. [7] an example of an extrapolation in which 
the plaquette decreases with A^, and in Fig. [8] an example in which it increases. Results are 
for = 10 — 53 and are plotted versus 1/A^. We use this variable because we have found that 
we cannot obtain a reasonable fit without taking the leading correction to be proportional to 
1/A^ (rather than 1/N'^). Examples of fits to Cq + Ci/N and Cq + Ci/N + C2/N'^ are shown in the 
figures, fitting either to all the data or dropping the two lowest values of A^. We find that fits 
of Co + Ci/N to the highest six values of A^ are tolerable (probability p ^ 0.04) for all choices 
b and k, that we have considered, and use such fits to obtain the results given in Table [H 

The results for the plaquette at A^ = oo are in striking agreement with the semi-quantitative 
prediction explained above. In particular, for k = 0.06 and 0.09 they are consistent with the 
pure-gauge large- volume results, while for n = 0.12 (close to Kc) they begin to differ. Thus, 
looking back at the scans of the plaquette in Figs. [3] and HI we see that the extrapolation to 
N = oo leads to an almost constant value between the onset of the funnel at ~ 0.05 and 
Kc, with the value being close to that of the pure-gauge theory. 

We also have results for a single point above Kc'- ^ = 0.15 at 6 = 1.0. We find here that the 
extrapolated plaquette differs, with high significance, from that below k^. This is consistent 
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FIG. 7. Plaquette vs. at b = 0.35, k = 0.09, showing various fits. 




FIG. 8. As in Fig. [7] but for hc = 0.12. 
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b 


K 


xVd.o.f. 


Cl 


Co 


pure-gauge value 


0.35 


0.06 


1.8 


0.75(4) 


0.549(1) 


0.550 


0.35 


0.09 


2.4 


0.34(4) 


0.552(1) 


0.550 


0.35 


0.12 


1.5 


-0.92(3) 


0.565(1) 


0.550 


1.0 


0.06 


0.2 


0.120(3) 


0.8694(1) 


0.8692 


1.0 


0.09 


1.1 


0.076(3) 


0.8697(1) 


0.8692 


1.0 


0.12 


0.6 


-0.248(4) 


0.8709(1) 


0.8692 


1.0 


0.15 


2.3 


0.39(1) 


0.8795(4) 


0.8692 



TABLE I. Results from large- extrapolation of plaquette expectation values. Extrapolations are 
done using a fit of cq + ci/N to results at N = 23, 30, 37, 42, 47 and 53. Results for cq, ci and 
X^/d.o.f. are given, with errors being statistical. Systematic errors (from different choices of fit 
function) are a few times larger than the statistical errors. Our best estimate of the pure-gauge 
large- volume expectation value is also quoted. The b = 1 value is obtained from Ref. [8|, while that 
at b = 0.35 is obtained from the N = 8 pure gauge result at 6 = 0.3504 from Ref. 



with our semi-quantitative prediction, and also indicates that the first-order transition at Kc 
survives the large limit. 

We now return to our numerical finding that the leading corrections to the plaquette scale 



as 1/A^. This result has also been found in the numerical studies of Ref. Uj. It differs from the 
naive expectation that, with fields in the adjoint irrep, corrections should be powers of l/N"^. 
We find, however, that fits to Cq + c^j^'^ + c^/N'^ are only possible with very large coefficients 
(|c2| ~ 10 — 20, \c/i\ ~ 2500) of alternating signs. We consider these fits unreasonable since we 
expect coefficients of 0(1). 

In fact, there are (at least) two possible sources of 1/A^ terms. The first can be seen from 



the perturbative result for the plaquette when one has a center-sjmametric vacuum jj] 

np = 1-1(1 -l/Ar) + 0(l/62), (4.7) 
which manifestly contains a 1/A^ correction]^ This correction arises from the fact that in the 



We thank Ari Hietanen for reminding us of this result. 
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decomposition fl4.5p non-trivial fluctuations in lie in SU{N)/U{1)^^^. In other words, the 
fluctuations must be off-diagonal, leading to the factor A^(A^ — 1) = iV^(l — 1/^)- We expect 
that the one- loop form (14. 7p should work reasonably well at 6 = 1 (as it does for cq in Table 
and that the predicted correction should be most applicable for the smallest values of k 
(where fermionic contributions to the plaquette are minimized). Indeed the result for ci at 
6 = 1, K = 0.06 lies close to the prediction of 1/8. 

A second source for corrections are the "would-be zero modes" of the Wilson-Dirac 
operator Dw, which we discuss in more detail in Sec. IV A[ There are 4(A^ — 1) of these (cor- 
responding, as in the gauge case above, to the diagonal generators of SU(N) in perturbation 
theory), and they form an 0{1/N) fraction of the total number of modes. Unless the con- 
tribution of these modes is exactly canceled by an 0{1/N) contribution from the remaining 
4(A^^ — A^) eigenvalues, these modes can cause observables to depend on odd powers of 1/A^. 
Our results for the spectrum of Dw suggest that they play an important role in the dynamics 
for the values of b at which we simulate. 

We now turn to the extrapolations of (|-P^P) and (|M^j,p). In the large- A^ limit, these can 
be written, using factorization, as |(-P^)P and \{M^^)\'^, respectively, both of which vanish if 
the center symmetry is unbroken. Thus an important test of our tentative phase diagram is 
that (l-P^P) and (|M^i,p) extrapolate to zero within the funnel. 

Ordinarily, corrections to factorization are proportional to l/N"^, but, in light of our expe- 
rience with the plaquette, we might also see 1/A^ corrections. In Fig. [9] we plot (|-Pip) versus 
l/N"^ for b = 0.35 and two values of k. There is qualitative agreement with a l/N"^ fall-off for 
both k's, but flts to a pure l/N"^ form, examples of which are shown in the flgure, have low 
confidence-levels. Satisfactory fits (one example of which is shown) can be found by adding a 
1/A^ term and dropping the lowest two values of A^. 

Results from such fits, for all the values of k and b at which we have done runs up to 
A^ = 53, are collected in Table ITTl The fits to (|-Pip) all have reasonable confidence levels. The 
coefficient of the 1/A^ term is small in all cases, and in fact is consistent with zero (within 
~ 3a) except for b = 0.35, n = 0.06. We also show results of fits to a constant plus 1/A^^ 
term. The fits are of very similar quality, and the constant turns out to be very small, and 
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FIG. 9. versus l/A^^ for h = 0.35 and k = 0.06 and 0.12, along with a variety of fits. 

consistent with zero except, again, at 6 = 0.35, k = 0.06. We conclude, aside from this one 
point near to the edge of the funnel, that the behavior of Polyakov loop is consistent with the 
hypothesis that reduction holds in the funnel. 

Turning to the corner variables, examples of which are shown in Fig. [10] with full results 
collected in the Table, we find a surprising result: at 6 = 0.35, k = 0.06, (|Mi2p) starts to 
increase once N exceeds 40, and clearly does not extrapolate to zero. The simplest interpre- 
tation of this result is that the center symmetry is broken for ^ 40. There is, however, no 
other evidence for such symmetry breaking. In particular, the distributions of the M^i, and 

are approximately uniform around the origin, the traces Kn of Eq. (14.21) are all consistent 
with zero, and the link eigenvalues (to be discussed below) are distributed uniformly. 

Instead, our favored interpretation is that the increase in (|Mi2p) with is due to the 
lower edge of the funnel, Kf, increasing with A^. This increase can be seen (albeit for 6 = 1) 
in Fig. [5] by comparing the results at = 16 and 30. It is quite possible that, for b = 0.35, 
as A^ increases, Kf approaches 0.06, possibly exceeding this value for A^ > 53. This would 
lead to the observed increase in (|Mi2p) since this quantity increases as one approaches the 
transition (as can be seen in Fig. [5]). This could also explain why our fits to (|-PiP) were less 
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0.35 


0.06 


0.006(1) 


1.21(6) 


0.36 


9(1) X 10-5 


1.31(1) 


0.39 


{\Pi\') 


0.35 


0.09 


0.0014(9) 


0.73(3) 


0.67 


2(1) X 10-^ 


0.76(2) 


0.73 


{\Pi?) 


0.35 


0.12 


0.001(3) 


0.57(3) 


0.82 


0(2) X 10-^ 


0.56(1) 


0.84 




0.35 


0.09 


0.152(5) 


2.9(1) 


0.42 


0.0023(3) 


5.3(4) 


2.57 




0.35 


0.12 


0.036(6) 


3.5(2) 


1.0 


5(1) X 10-"* 


4.1(2) 


1.2 


{\Pi\') 


1.0 


0.06 


-0.0001(3) 


1.17(1) 


0.025 


-1(3) X 10-6 


1.17(1) 


0.025 


m') 


1.0 


0.09 


-0.0003(4) 


0.70(1) 


1.2 


-5(6) X 10-6 


0.70(1) 


1.2 


m') 


1.0 


0.12 


-0.0010(3) 


0.55(1) 


0.60 


-1.4(4) X 10-5 


0.54(1) 


0.58 




1.0 


0.06 


0.69(3) 


-1.6(7) 


0.58 


0.010(1) 


9(2) 


2.7 


{\Ml2\') 


1.0 


0.09 


0.0053(7) 


6.1(2) 


0.99 


8(2) X 10-'' 


6.9(3) 


1.7 




1.0 


0.12 


0.01(1) 


5.5(4) 


0.4 


0(1) X 10-* 


5.4(2) 


0.41 



TABLE II. Results from fits to the large-iV behavior of {\Pi\^) and (|Mi2p). Fits are to = 23, 30, 
37, 42, 47 and 53 using fi{N) = ci/N + C2/N'^ and /2(iV) = Cq + ^/A^, and the quoted errors are 



statistical. 

satisfactory at b = 0.35, k = 0.06. 

For all the other values of k that we have considered (|Mi2p) decreases monotonically with 
A^. This is exemplified by the k = 0.12 results in Fig. [TOl As for the Polyakov loops, a pure 
1/A^^ fit fails in most cases, but here we find (cf. Table [IT|) that the addition of a 1/A^ term 
usually leads to a better fit than the inclusion of a constant. We also find that, in almost 
all cases, the required 1/A^ (or constant) term has a coefficient which differs significantly 
from zero. We have also done fits to C2/N'^ + c^/N^^ (an example is shown if Fig. [TOl) but 
the fits require very large coefficients having opposite signs, a fine-tuning which we consider 
unlikely to be the correct description. Overall, we think the most reasonable fits are those to 
Ci/N + C2/A^^, because they have the highest confidence levels, and because we have seen in 
the plaquette that 1/A^ terms are needed. 

The results presented so far are consistent with the funnel (in which center symmetry 
is unbroken) remaining of finite width as A^ — )■ 00, so that reduction holds for masses up 
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FIG. 10. As for Fig. [9]but for the corner variable |Mi2p. Only fits to k = 0.12 are shown. 

to 0(l/a). We cannot definitively draw this conclusion, however, because of the following 
scenario. Imagine that the funnel width vanishes (for any fixed 6) as — ?■ oo. Then, for each 
b, K point in the putative funnel, symmetry breaking would occur at a (possibly large, but) 
finite, value of A^. Nevertheless, there would be a finite range of for which the symmetry is 
unbroken, within which the arguments for reduction hold. Appropriate variables (such as the 
plaquette) would equal infinite volume values up to corrections proportional to powers of 1/A^. 
Thus, within this range, it might appear that one can extrapolate to the symmetry-unbroken 
N = oo limit, but this would in fact not be the case. The results for the plaquette at 6 = 0.35, 
K = 0.06 in Fig. IHlare an example of such misleading scaling, since we have strong evidence 
from the corner variables that Hf > 0.06 at this b for large enough A^. 

D. A^-scaling of the funnel width 

In the light of the results in the previous subsections it is important to directly study the 
A^ dependence of the funnel width. We have done so by focusing on Kf, the lower edge of the 
funnel. If we can show that Hf remains below as A^ — > oo, then the funnel remains open. 
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FIG. 11. Fine scans of absolute values of all Polyakov loops and corner variables in the low k region, 
for AT = 30 at 6 = 0.75. 

To investigate this issue we have done fine scans of the small k region, an example of 
which is shown in Fig. [Ill There are two phases before one enters the funnel: a Zi phase for 
< K ^ 0.02 and a Z2 phase from 0.02 ^ n ^ 0.05. The transition between the first and 
second phase shows significant hysteresis, while that between the second phase and the funnel 
does not. 

Determining K,f{N) to high precision is a significant numerical challenge. We have thus 
focused on a single value of coupling, 6 = 1. For this b, Fig. [3] shows that hes in the range 
0.13 — 0.14. We have done very fine scans near the edge of the funnel (roughly n = 0.02 — 0.07) 
with up to 53, and find that the corner variables are the most useful in determining the 
transition. We are able to pin down the transition, conservatively, to about 6k = ±0.001. The 
transition from the funnel is to a Z2 phase for most (as in Fig. [TTl) . but to a Z3 phase for 
N = 47 and 53. The results are plotted against 1/A^ in Fig. [12], and show remarkable linearity 
(note that, as earlier, we have excluded = 10 and 16 — including them requires adding a 
l/N"^ term to obtain a satisfactory fit). Two fits are shown. The first is to Cq + Ci/N, and has 
a very small x^/d.o.f., with a reasonable coefficient of 1/A^. It yields K,f{N = 00) = 0.0655(5), 
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FIG. 12. The dependence of (lower edge of the funnel) on for 6 = 1. The fit functions are 
discussed in the text. 

a value far below Kc- The second fit is to Kc + Ci/N + C2/N^, with the intercept fixed to 
He = 0.125. This fit is extremely poor, and it gets even worse for Kc = 0.13 — 0.14. We 
have also done the corresponding fits to the alternative quantity aruf = 1/2k/ — 1/2kc, using 
Kc = 0.125 or 0.13, and find consistent results. We conclude that^at least at this value of b, 
the funnel has finite width when N = oo, so that reduction holdscl 

E. Distributions of link eigenvalues 

We find that histograms of link eigenvalues provide very useful information on the nature of 
symmetry breaking on either side of the funnel. They are also sensitive to patterns of symmetry 
breaking in which both Polyakov loops and corner variables vanish, and thus provide a more 
stringent test that the symmetry is indeed unbroken in the funnel. In this section we present 
examples of the results that allow us to fill in the details of the phase diagram of Fig. [2J 
Note that, if we use the linear fit, then the funnel at 6 = 1 passes k — 0.06 when N « 92. Thus the successful 
extrapolations of the plaquette, and (|Mi2p) for = 1, k = 0.06, presented in Tables HI and HIl are 

examples of the phenomenon described above in which reduction only holds for a window of values of N. 
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We begin with an example of a histogram within the funnel, shown in Fig. 13(a) The 
eigenvalues are taken to have the range — vr < < vr, and are collected in 3A^ bins of width 
27r/(3A^). Thus symmetry implies invariance under periodic translations by multiples of 
3 bins. In fact the distribution is consistent with being uniform]^ 

When we move outside the funnel the attraction between eigenvalues leads to formation of 
clumps in the complex plane. The number of clumps k identifies the approximate remnant 
symmetry and generally decreases as we move away from the funnel. In Fig. [12] we provide 



several examples of clumping patterns. Figs. 13(b) and 13(c) present Zi and Z2 phases in 



the small k region while Figs. 13(d), 13(e) and 13(f) show Z3, Z4 and Z^ phases in the large 
K region. Note that the partial symmetry breaking can also be seen in the corner variables 
giving complex patterns (compare Fig. [5]). Polyakov loops are much less sensitive to this 
partial symmetry breaking, since they almost vanish due to the approximate Z^ symmetry. 

we 



We find that the remnant symmetry is not always exact. For example, in Fig. 13(c' 



have two clumps for = 23 and in Fig. 13(d) we have three clumps for N = 16. Therefore 
the eigenvalues cannot be equally distributed between the clumps and the symmetry is only 



approximate. Even in Fig. 13(f) which shows five clumps for = 30, we see that the clumps 
are not even and correspond to 7,6,7,4,6 eigenvalues, respectively. We also find that different 
runs can have different patterns of eigenvalue clumping, e.g. 7,7,6,6,4 versus 7,6,6,6,5, but that 
it is rare for the clumping to change during a run. Thus it appears that there are competing 
"vacua" which are not exactly related by center symmetry transformations. 

To fully understand the symmetry breaking, we need to know whether there are correlations 
between the eigenvalues of different links. What we find is that, whenever the center symmetry 
is broken, the eigenvalues for all four links are highly correlated. To illustrate this, we use 
a case with three clumps which makes the results easy to visualize. Figure [TH shows the 
resultant clumping and correlations. Here we apply a gauge transformation which diagonalizes 
Ui and orders the phases 6f, and then look at the phases of the diagonal elements of f/2,3,4- 
These matrices are close to diagonal, so these phases are presumably close to those of their 
eigenvalues. Recall that there is no ambiguity in the ordering of the diagonal elements once 
As a check on our code, we have calculated the distribution in the b — 0, k — limit, i.e. for the Haar 
measure on the links, and obtain the theoretically expected form, which is ZAr-invariant, but does oscillate 



within each segment, although the amplitude of the oscillations falls as 1/A'' 
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(b) N = 23,b= 1.0, K = 0.01, 
2000 configs, all links 



(d) iV = 16, 5= 1.0, K = 0.24, 
5000 configs, C/3 only 




(f) TV = 30, 6 = 1.0, K = 0.23, 
1000 configs, Ui only 



FIG. 13. Histograms of the phases 9^ of the link eigenvalues. More details of the binning are 
discussed in the text. 
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FIG. 14. Results for the phases of the diagonal elements of for 20 thermalized configurations at 
N = 16, b = 0.35 and k = 0.23 on an UP scan. Phases are determined in a gauge such that Ui is 
diagonal with the phases ordered. For further discussion, see text. 



we specify the order for Ui. The result shows that the clumps (of 6, 4 and 6), while being 
positioned at different angles, are almost completely correlated between all four links, and do 
not change during the Monte-Carlo evolution. Because of these correlations, the approximate 



remnant of the center symmetry for the parameters of Figs. 13(d) and[T3]is Z^, and not Z|. 

Once outside the funnel, the number of clumps decreases as we move to higher k. We have 
extended some runs to k = 0.6 and find that the UP scans end up in a two clump state, while 
the DOWN runs, which begin from an ordered start, begin with a single clump, and then have 
a transition, as k, is decreased, to two (well separated) clumps. In fact, the transition appears 
to occur in stages where more and more eigenvalues peel off from the initial clump. As k is 
further decreased, the number of clumps increases until we enter the funnel and there is no 
longer any clumping. The largest number of clumps depends on A^, and the largest we have 



observed is five, as shown in Fig. 13(f) 



The changes in clumping for the large k values, surveyed above, are qualitatively consistent 

28 



with the arguments presented in Ref. ll|] based on the one-loop free-energy for the hnk eigen- 
values. Decreasing k from a large value corresponds to reducing the quark mass am. For large 
am, gluonic interactions dominate the free energy, and lead to attraction, and thus a single 
clump. As am is reduced, the fermionic contributions lead to repulsion at large eigenvalue sep- 
arations (corresponding to large momenta, so that the mass term is unimportant), while there 
remains attraction for small separations. This allows the possibility of two clumps. Reducing 
am further the repulsion becomes important for smaller eigenvalue separations, and clumps 
are pulled apart into a greater number of stable clumps. At the same time, quantum fluctua- 
tions within each clump are always present, so that the clumps have a finite width (which is 
proportional to b^^^^ for weak coupling). Eventually, as the number of clumps increases, the 
distance between the clumps is smaller than the widths, and the clumping is washed out. 

We would expect a similar sequence of dumpings to occur as we increase k from zero, since 
this also corresponds to reducing am. This is indeed what we observe, although the maximal 
number of clumps is smaller on this side of the funnel. For < 23 we only find a Zi phase, 
for 23 < A^ < 47 we see both a Zi and Z2 phase (see Fig. [TTl) . while for A^ = 47 and 53 we 
find Zi, Z2 and Z^ phases. At larger values of b, the arguments of Ref. ll| imply that the 
maximal number of clumps should increase. Indeed, we do find that, as b increases, the Z3 
phase appears at smaller values of A^. 



F. Results at large b 

The perturbative calculations of Refs. jol, Q lead us to expect that the center-symmetry- 
unbroken funnel will close as 6 —t- 00, so that, in the continuum limit, reduction only holds 
for m = (for A^ 00). In addition, Ref. [lli makes a prediction for how rapidly the 
funnel should close: its width in am (and thus in k) should be proportional to b~^^^. This is 
because the width of each clump of link eigenvalues is predicted to scale to zero proportional 
to 6^^/^. A second prediction (which is explained in the previous subsection) is that, at the 
edge of the funnel, there should be multiple phases with differing numbers of clumps, and 
that the maximum number of clumps should increase with b (as long as A^ is large enough). 
In other words, a behavior similar to that we have already seen on the right side of the funnel 
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FIG. 15. The average plaquette in scans at extremely high b for N = 10. Note the highly compressed 
vertical scale at large b. 



{Zi — i- ^2 —i- ■■■—)■ Z5; see Fig. [2]) extends to larger groups as b increases. 

We have investigated these predictions by doing scans in n for = 10 (and, in some cases 
= 30) at 6 = 5, 10, 50 and 200. We find that the HMC algorithm mostly performs well even 
at these very weak couplings. We did have to reduce the step-size as b increases, such that 
^A/z)/acceptance increases roughly like y/b. We also find that, for b = 200, thermalization 
for each new value of k sometimes takes longer than our allotted 450 trajectories. On the 
other hand, the number of CG iterations gradually decreases with increasing b. We also note 
that run histories of observables show no indication of correlation times that are close to the 
number of trajectories we were using for measurements (7500). 
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FIG. 16. The dependence of Hf on b for = 10 and 30. The vertical (blue) line marks 
Kc{b = oo) = 1/8. The fit functions are discussed in Sec. IIVFI 

Results for the plaquette are collected in Fig. [151 The general shape of each curve is similar 
to those at 6 = 1 (see Fig. [3]) but the jump at the putative Kc falls rapidly with increasing b. 
This is qualitatively consistent with the expectations from chiral perturbation theory if this 
is the first-order scenario of Ref . l20|l . 

To venfy the Kvpo.hes. of Ref. Q *Kat an., ^ 1/2.. - 1/2. oc „e a.a.y.e the 
lower edge of the funnel as a function of b. It would obviously be advantageous to repeat the 
analysis of Sec. IIVDI for all values of 6; unfortunately this is numerically too demanding. We 
do, however, have estimates of at = 10 and 30 for a wide range of b. These are shown in 
Fig. [161 We find that it is harder to determine k j for b away from unity. For smaller values, 
e.g. b ^ 0.35, the transition becomes smoother. For much larger values, there is significant 
hysteresis (as seen in the middle panel of Fig. [T5l) . The net result is that the errors in are 
much larger than those at 6 = 1. 

Figure [THl also shows fits to anif = c b~^^^, with Kc set to 0.125 for all values of b for 
simplicity. The fit at = 10 is good, while that at = 30 is poorer. Better fits at A^ = 30 
can be obtained using estimates for the actual value of at each b, but these estimates have 
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sufficient uncertainty that the resulting errors in mj are substantially increased, so that the 
agreement with the theoretical form is less significant. Overall we conclude that our results 
are consistent with the predicted dependence on b. 

The figure also indicates that the narrowing of the funnel as increases holds for all values 
of b. We do not attempt to extrapolate to = oo from = 10 and 30, as our experience 
with 6=1 indicates that = 10 is not in the asymptotic 1/A^ region. However, given that 
we do find a finite width when N ^ oo aX b = 1, the observation of mild dependence on b 
suggests that the funnel will remain of finite width also at other values of b. 

V. MEASUREMENTS INSIDE THE "FUNNEL" 

In this section we make a detailed study of the funnel region in which reduction appears to 
hold using results from A^ = 10 — 53. We consider in turn the spectrum of the single-site Wilson 
operator Dw, the spectrum of Q^, and, finally, attempt to extract a physical observable — the 
heavy-quark potential — from large Wilson loops. 

A. Spectrum of Dyy 

One way of viewing the equivalence of single-site and large-volume theories is that the 
space-time volume is being packaged inside the gauge matrices. It is thus useful to introduce 
an effective size. Lea, and corresponding effective volume, Lg^, and study their scaling with 
A^. What we mean by L^s is that the single-site theory leads to the same physical results 
as the theory on an volume with N^s colors, up to corrections suppressed by powers of 
1/Nes- Clearly there is a trade-off between increasing L^s and increasing A"eff. Here we take 
the approach of holding A^efr fixed, but large enough that 1/Ncs corrections to quantities of 
interest are small, and then asking how Lcs scales with A^. 

Within this framework, the most conservative possibility is provided by orbifold-based 
demonstration of volume independence This demonstration also provides an explicit 
example of the packaging of the volume into the gauge matrices: the N x N link matrices 
are partitioned into blocks of size N^s x A^es, with A'efi = N/L^^, only L^g of which are 
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non-zero, and the resultant orbifolded theory is an SU {N^s) gauge theory on an L^g lattice. 
Equivalence is demonstrated by holding Les fixed, and taking N, and thus also A'efr, to infinity. 
Instead, using the approach espoused above, if we hold N^a fixed, then increasing N leads to 



A less conservative possibility is obtained, following Refs. [33|, |3J], by assuming that all 
entries in the link matrices are used in the packaging of the large volume theory (not just 1 
out of every L^q as in the orbifold construction). This leads to L^g ~ A^^ or L^s oc N^^"^. 
There is also a more optimistic possibility, L^s oc A^, which is motivated in the Appendix. 

The spectrum of the fermion matrix, D^, can help distinguish these possibilities, as well 
as give insight into the nature of corrections to the large- A^ limit. We expect, if reduction 
holds, that the spectrum should resemble that of a large-volume four-dimensional theory on 
an Lgg lattice. In particular, for weak couplings, 6^1, the spectrum should have the familiar 
five "fingers" which reach down to the real axis. The number of fingers is a direct indicator 
of the dimensionality (there are d+1 in d dimensions), and the distance of the eigenvalues in 
the fingers from the real axis should scale as l/Leg. These points are discussed in more detail 
in the Appendix. 

We now show some representative results for the spectrum of Dw^mQ = 0) from our 
simulations. The operator in the determinant is 

-^[4Dw{0)+mo] , (5.1) 
4 + mo 

so that eigenvalues of 4Dvy(0) close to A = —rriQ = 4 — 1/2/t are suppressed. Since the spectrum 
is bounded, < ReA < 8, the determinant suppression is important only for n > 1/8. We 
also note that, unlike on a lattice with an even number of sites in each direction, the spectrum 
is not symmetric under refiection about the ReA = 4 axis. Thus the first (the leftmost) and 
fifth fingers are not related by symmetry, and neither are the second and fourth. If such a 
symmetry holds approximately, it indicates the presence of reduction. 

In Fig. [T7]we show how the spectrum changes as we vary k, at fixed b = 1.0 and A^ = 16. 
At K, = 0.01, where we are in the Zi phase (see Figs. ElandE]), we see one main finger and a 
small indication of a second. This is consistent with the eigenvalues forming a single clump, so 
that the momenta, given by eigenvalue differences, are all small. At k = 0.03 we have moved 
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into the Z2 phase, with two clumps of eigenvalues. We see that this allows the spectrum 
to spread into all five fingers, because eigenvalue differences can now range up to tt. The 
details of the spectrum differ from those of a large- volume free fermion, however, in particular 
having a low density of points in the central three fingers and a "rectangular" shaped envelope. 
Nevertheless, it is clear that one must interpret the spectrum with care — the presence of five 
fingers alone does not imply that reduction holds. 

The next value, k = 0.12, is well inside the funnel, and we see a distribution which is 
qualitatively similar to that of a free fermion, with a rounded top and five fingers. These 
features are present for all k < Kc inside the funnel. Particularly noteworthy is the presence 
of the comet-shaped clump of eigenvalues near the origin. We find that there are exactly 
4(A^ — 1) eigenvalues per configuration in this clump. We thus interpret them as would-be 
zero modes, i.e. eigenvalues that would be zero if 6 — > 00. These modes are dropped in weak 
coupling calculations, both because they do not impact the dynamics (as they do not depend 
on the and because they form only an 0{1/N) fraction of the total number of modes. The 
spectrum indicates, however, that they could have an important impact on the long-distance 
dynamics which might overcome their relative paucity. We recall that for a large-volume 
Wilson operator it is the small eigenvalues which determine long-distance behavior such as 
chiral symmetry breaking. For very large N we expect small eigenvalues to come dominantly 
from the first finger, which should approach the real axis. What we see from the figure is that 

= 16 is quite far from this limit. Thus we conclude that the would-be zero modes are a 
potential source of the 1/N corrections observed above in the plaquette and other quantities, 
and that their contribution could be sizable (given how far the "true" low-energy modes in 
the first finger are from the real axis) . 

The spectrum within the funnel but just above the transition is illustrated by the result 
for K = 0.14. Eigenvalues near ReA = 0.43 are suppressed by the determinant. The would-be 
zero modes cluster to the left of this excluded point, while the first finger now approaches 
closer to the real axis. The latter feature indicates that the funnel-region above the transition 
is more continuum-like, which is consistent with its larger average plaquette. On the other 
hand, the spectrum as a whole is less symmetric about the ReA = 4 line than that below the 
transition. 



34 



3 4 5 6 7 

Re I 



4 
Re?i, 



(a) K = 0.01 



(b) K = 0.03 




3 4 5 6 7 

Re A 



(c) K = 0.12 





(d) K ^ 0.14 




(e) K = 0.17 (f) K = 0.24 

FIG. 17. Spectrum of 4Z)p^/(mo) from simulations at 6 = 1.0 and = 30 at k = 0.01, 0.03, 0.12, 



0.14, 0.17 and 0.24. Only eigenvalues with positive imaginary part are shown. 
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FIG. 18. Spectrum of 4Z)vi^(mo) at 6 = 1.0 and k = 0.12, for = 37 and 53. Note that since 
(53/37)^ ~ 2 the number of points is approximately the same in both plots. 

Moving to K = 0.17, which is still inside the funnel, the would-be zero modes have spread 
out again (perhaps because the excluded point has now moved to A = 1.06), while the first 
finger has become longer and denser. The second finger, however, has almost disappeared. 

Finally, at k = 0.24 we are in the phase. This is reflected by the spectrum breaking 
into three distinct regions (only two being visible since the third has negative imaginary part), 
resulting from eigenvalue differences distributed around and ±27r/3. 

We have done similar scans at lower b, but the results are less illuminating, because the 
bulk of the spectrum moves closer to the real axis, such that, at 6 = 0.35, one cannot see 
any fingers. A better approach is to use the spectrum of Q^, as described in the following 
subsection. 

We have also studied the N dependence of the spectrum at 6 = 1 and k = 0.12. Results 
from = 37 and 53 are shown in Fig. [181 ^ind can be compared to the = 16 results in 



Fig. 17(c) The spectra at A^ = 37 and 53 differ very little. The main changes are that the 
size of the clump of would-be zero-modes decrease as A^ increases, and that the tip of the 
first and fifth fingers move down slightly. The tips of the other fingers, however, barely move. 
Compared to A^ = 16, on the other hand, the fingers are somewhat more extended. 

We draw several conclusions from these results. First, the qualitative agreement of the 
spectrum within the funnel with that from a large volume four- dimensional theory supports 
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our conclusion that reduction holds therein. Second, the Lcs ~ A^^^^, crystalline distribution 
of eigenvalues described in the Appendix is disfavored, since, for our values of N, it would 
not lead to the presence of fingers, and thus differ from our resultso Third, our results 
are also inconsistent with the Leg ~ N model, since the fingertips do not approach the real 
axis as fast as the expected Fourth, the would-be zero- modes are a possible source of 

0(1/A^) corrections. And, finally, these zero-modes may make an important contribution to 
dynamics, thus enhancing the 1/A^ corrections (again, provided that these 1/A^ corrections 
are not exactly canceled by 1/A^ corrections from the 4(A^^ — A^) modes in the bulk of the 
eigenvalue distribution) . 

We have also calculated the spectrum of the Dirac operator in the fundamental representa- 
tion. This gives information directly about the link eigenvalues, rather than their differences. 
The results confirm our understanding of the phase diagram and eigenvalue distributions 
explained above, but are not shown for the sake of brevity. 



B. Spectrum of 

An alternative view is provided by the spectrum of the squared hermitian Wilson-Dirac 
operator, = _Dvi/(mo)-DvK(mo)^ Its eigenvalues, Xq2, are real and positive. In the contin- 
uum limit, the spectrum has a gap, turning on at Aq2 = (arriphys)^. Away from the continuum 
limit, the turn-on is smoothed, but still begins approximately at the square of physical bare 
quark mass jssj . For small enough quark masses, and if there is spontaneous chiral symmetry 
breaking, the spectrum above the gap is approximately constant, with a value proportional to 
the condensate. Thus the spectrum can teach us about the size of the quark mass and about 
long-distance physics (assuming reduction holds). The information is also contained in the 
spectrum of Dw{0), but for 6 < 1, when the fingers are obscured, is hard to extract. Thus we 
have used the spectrum of mainly for b = 0.35, which, we recall, is a bare coupling close 
to those used in typical large-volume simulations. 

Results for b = 0.35, k = 0.12 and A^ = 10 — 47 are shown in Fig. [T9l (Results at A^ = 53 are 
We have studied this further by calculating the spectrum on a 2^ lattice with iV = 3, corresponding to a 
"partial crystallization" of a single-site iV = 48 theory. This spectrum also has no fingers. 
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FIG. 19. Spectrum of Q"^ aX. h = 0.35 and k = 0.12, for = 10, 20, 37 and 47, using 150, 150, 20 
and 60 configurations, respectively. The vertical scale is arbitrary, but the relative scales for different 
N are chosen so that the area under each spectrum is the same. Errors are not shown, but can be 
estimated from the kinks in the spectra. 

very similar to those at = 47 but have lower statistics and are thus not shown.) The spectra 
are normalized to have the same integral, so that the large limit can be taken. The peak 
at small eigenvalues has the correct area to contain just the would-be zero modes, so that its 
area drops as in the normalized spectrum. It will disappear entirely when — t- oo. The 
bulk of the eigenvalues form a "hump" which at the upper end (Aq2 ^ 1.5) is approximately 
independent of A^, while at the lower end depends on A^. The form of the hump at A^ = 47 
should give a good approximation to the form at A^ = oo because the small-eigenvalue peak 
has httle area left to "redistribute" to the hump. A crude extrapolation of the leading edge of 
the hump at A^ = 47 suggests that the gap at A^ = oo will be at Aq2 0.1. This corresponds 
to a quark mass of amphys ~ 0.3, modulo the unknown renormalization factor, which, however, 
we expect to be of 0(1). This is the only "measurement" of the quark mass that we have 
obtained, and shows that the quark is relatively heavy, not much below mphys = 1/a. 
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These results shed hght on the issue of how large N needs to be for reduction to be useful. 
On the one hand, one can seen that the bulk of the eigenvalues (those in the hump) are close 
to their large N limit by ~ 40. On the other hand, the would-be zero modes, though 
making up only an 0{1/N) fraction of the total, are the dominant contribution in the low 
mode region. Which effect is more important is not clear a priori — one must calculate physical 
observables and study their iV- dependence. 



C. Large Wilson loops 

Our ultimate aim in studying the single-site models is to use them to calculate physical 
quantities in phases where large-N reduction holds. An important quantity that should be 
accessible in such phases is the heavy-quark potential. To obtain this we calculate rectangular 
Wilson loops using the large-N reduction recipe (and also averaging over orientations): 

WiLuL,) = ^Y.(^^^^' {U;:^U^'Ul'^^Ul'^^)) , (5.2) 



12 ^\N 

For N oo the result should equal the infinite-volume large- value. The potential can be 
obtained as usual from the large behavior 

(Li, L2) c(Li)e-^(^^)^^ . (5.3) 

At large Li we expect linear behavior if we are in a confining regime: 

^^^-a. (5.4) 

For finite A^, reduction will only give useful results for loops whose sizes satisfy Lj <^ N, 
and the key question is how much smaller than A^ do the Lj need to be. Another important 
question is whether it suffices to calculate Wilson loops using unsmeared links, i.e. whether 
the statistical errors will overwhelm any signal of interest. State-of-the-art calculations in 
large volumes use smearing, as well as other noise reduction techniques. 



In our earlier study of the Nf = 1 model js], in which A^ < 13, we did not find a useful 
"window" as a function of the Lj. For fixed Li, for example, the dependence on L2 was a rapid 
drop for a few points followed by a slow rise. The drop was not extensive enough to determine 
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the potential from the coefficient of the exponential Here we have results extending up to 

= 53, and thus expect that the situation will be substantially better. 

In Fig. [20] we show results for 1 x L loops for a range of values of on a log-linear scale. 
For each value of we find an approximately exponential decrease followed by slow, roughly 
linear increase. The latter we interpret as a finite effect, since it begins at larger values of 
L a.s N increases. The good news from this plot is that we see convergence to the expected 
exponential drop-off as A^ increases. For example, at L = 6 the A^ = 37 point has peeled off 
the linear envelope, but the N = 47 and 53 points are in good agreement. This is large-A^ 
reduction in action. The bad news is that the maximum value of L at which convergence 
occurs, Ljnax, giows ouly slowly with A^. This is not unexpected: we are trying to extract an 
exponentially falling contribution to a quantity which has finite A^ corrections. To estimate 
the size of these corrections, one can look at the the minimum value of the loop as a function 
of L, which we find falls as approximately 1/A^. Thus Lmax grows only logarithmically with 
A^, which presents a significant numerical challenge. 

Despite this challenge, we see from the figure that we can extract a value for the slope 
at small L with reasonably small errors. This gives —1^(1), the potential at unit separation. 
To extract a, we need the potential at larger separations. We show in Fig. EH the results for 
5 X L loops. The overall pattern is similar to those for 1 x L loops, but the convergence in 
A^ of the falling parts of the curves is much poorer. Only results for L < 2 appear converged. 
Thus we cannot extract ^(5) with any reliability. It is of course not a surprise that difficulty 
of determining V{L) increases with L, since the signal falls off more quickly while the 1/A^ 
background is little changed. 

Concerning statistical errors, we see from both of these plots that 20 configurations is 
sufficient to pull out the rapidly falling part of the curves. The problem is not the statistical 
errors, but rather the 1/A^ corrections. 

We have carried out a similar investigation at other points in the funnel. We find that as b 
increases, the slope of the initial fall-off decreases. If we move to the other side of the transition 
(where the average plaquette is closer to unity) we find that the slope decreases further, and 
also that the large-L approximately linear rise changes to an almost L-independent behavior. 

It would clearly be of interest to gain some understanding of the large L behavior of the 
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FIG. 21. As for Fig. [20l except for 5 x L loops. 
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Wilson loops. The only analytic approach that we are aware of that can shed some light on 
the issue is to calculate the loops in strong coupling. In the 6 = limit, and ignoring the 
fermion determinant, the gauge links are distributed according to the Haar measure, and one 
can show thali^ 

WiL,,L2)'^' j;^^{L, + L2-l-L,L2/N^) for < L^, < N . (5.5) 

For Lj ~ 0(1), the contribution is of 0(1/A^^), rising to of 0{1/N) when at least one of the 
Lj is of 0{N). This shows how the correlations can lead to an increasing result as one of 
the Lj is increased, which is qualitatively in agreement with the results of Figs. [20] and EH 
Quantitatively, this model does not, however, reproduce our data at 6 ~ 1. For example, 
we would expect that Wmm ~ l/N'^ in the model, but we find instead an approximate 
dependence. 



VI. DISCUSSION AND OUTLOOK 



We have presented a detailed study of the single-site version of large- QCD with two fla- 
vors of adjoint Dirac fermions — the so-called Adjoint Eguchi-Kawai (AEK) model — discretized 
using Wilson's gauge and fermion actions. This seemingly simple model turns out to have 
a rich phase structure, as shown in Fig. |2l some aspects of which can be understood semi- 



quant it at ively 111]. Our most important result is that we find, for up to 53, and for b 
up to 200, a broad funnel in the (/?, b) plane in which the Zff center symmetry is unbroken. 
This region encompasses both light {k — )■ Kc) and heavy (m ~ 1/a) quarks. While the funnel 
narrows as — ?■ oo, we present strong evidence that it remains of finite width in this limit. 
If so, then the single-site theory, when simulated within the funnel, is equivalent to the corre- 
sponding large-volume theory, up to corrections, suppressed by powers of 1/A^, which can be 
made arbitrarily small. Thus one can use the single-site model to study the large- A^ version 
of the minimal walking technicolor (MWT) model discussed in the Introduction. 

In particular, the phase structure we find within the funnel should be identical to that 
of the large-A^ MWT model. We find a first-order transition, with a discontinuity in the 

Note that these loops ah vanish in infinite vohime at strong coupHng. One obtains a non-zero result on a 

single-site lattice because the C/^ and C/^ terms in the loop are correlated. 
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plaquette (and other variables), for values of b at least up to 6 ~ 1 (see Fig. [T^ . This is in 
contrast to the results of direct large-volume simulations of the N = 2 two-adjoint model, 
which find that the transition changes from first-order at strong coupling to second-order at 



weak coupling, with the transition occurring at 6 0.25 |25l . |26| . The quark mass vanishes 
along the second-order portion of the transition line, and it is in part by studying the non- 
perturbative /3-function along this segment of the line that evidence has been found for an 
infrared fixed point. Further evidence that the N = 2 theory is conformal in the infrared 
comes from studies of the spectrum and other quantities as one approaches the second-order 
line (as summarized in Refs. 



15 



16|). 



The most straightforward interpretation of our phase diagram is that the theory is confining 
in the infrared, with chiral symmetry spontaneously broken, and that the first-order transition 



is an example of the first-order scenario of Ref. |20| in which the chiral condensate flips sign. 
The transition is a result of competition in the effective potential for the condensate between 
terms proportional to the fermion mass m and discretization errors of size a^A^. Here A is the 
confinement scale, and it is essential for this picture that such a scale is present. The anal- 
ysis predicts that the discontinuity should drop rapidly as b increases, which is qualitatively 
consistent with our numerical results. 

An alternative interpretation is that the transition is a "bulk" transition, which happens 
to extend to large b, but at some finite be > I becomes a second-order line. This would be an 
extreme version of what happens for the N = 2 theory, where the bulk transition extends to 
b ~ 0.25. In this picture discretization errors do not allow one to approach the massless theory, 
except for extremely weak coupling, so that one cannot determine the infrared properties of 
the massless theory!^ We view this as an unlikely possibility, since bulk transitions occur 
generically at strong coupling. 

Whatever the interpretation, it is striking that there is such a large difference in the phase 
diagrams of the N = 2 and N = oo theories. On the one hand, such a difference is in 



conflict with the models of Ref. 



for the position of the conformal window, which give 



N-independent predictions. On the other, we note that the beta-function (for the coupling b) 
We note that there is no problem in principle with using single-site models to study theories in the conformal 
window, as stressed in Ref. [12 1. 
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does depends on A^, starting at four-loop order. The situation clearly warrants further study. 

Another result of our work is that we find strong evidence for corrections when ex- 
trapolating the average plaquette to = oo. This is consistent with perturbation theory, 
which predicts effects at one- loop order for a single site theory, in contrast to the l/N"^ 
corrections that one finds in infinite volume. Based on our study of the spectrum of Dw-, the 
Wilson-Dirac operator, we suggest that another source of corrections may be the contri- 
butions of the 4(A^— 1) modes of Dw that become zero-modes when h — > oo. One caveat with 
our suggestion for the source of the behavior is that the contribution of the would-be 
zero-modes can be canceled by those from other modes. This happens in infinite-volume per- 
turbation theory (as in the results for the first two terms of the /3-function mentioned above) 
and also at strong coupling [as in the result for Wilson-loops of fixed size at strong coupling, 
eq. (15. 2p ]. We also note that the presence of significant dependence in the plaquette may 
be related to the apparent difference between the phase diagrams of the N = 2 and N = oo 
theories. 

If, as our results suggest, reduction holds within the funnel, then a lattice theory with 
multiple sites is being packaged inside the four link matrices and single-site fermion fields. 
One can think, approximately, of an effective lattice size. Leg. It is important for practical 
applications to determine how L^s grows with A^. Our results for the spectrum of (which, 
within the funnel, looks qualitatively similar to that on a large volume, cf . Fig. [T8|) suggest that 
neither the most pessimistic {Les oc N^^^) nor the most optimistic {Lcs oc A^) possibilities hold. 
Let us assume, then, that L^s oc N^^'^ — an intermediate possibility motivated in Sec. IV Al — 
and consider the question of whether, if one wants to determine the large- A^ properties of a 
theory, it is computationally advantageous to use a single-site theory or one on a large volume, 
L"^. In the former case, we have found from our simulations that, for fixed b and k, and with 
K near to Kc, 

CPU(l-site) oc A^^-^ oc L^f^N^ , (6.1) 

where in the final expression we have used L^s oc N^^'^. For the large volume theory we expect 
(for fixed lattice spacing and fermion mass) 



CPU(L^) oc L^A^^^ , 
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(6.2) 



where the is the standard hybrid Monte-Carlo volume scaling 37|, while N'^ is the operation 
count for the core operation of N x N matrix multiplies. This comparison suggests that the 
single-site approach could be computationally advantageous. While there are many caveats 
to this conclusion (e.g. the single-site scaling form is based on simulations for a finite range 
of N, and the scaling of Les with is not established), we take it as motivation to further 
pursue studies of reduced models. 

One aspect of such studies is calculating physical quantities such as the string tension and 
particle masses. We have taken a first step in this direction by calculating large Wilson loops 
and attempting to extract the heavy-quark potential. We find that we can do so for small 
separations (roughly out to 3 lattice spacings for b = 0.35) but not beyond. The difficulty 
arises because the signal must be determined from an exponential decay as a function of 
the loop size, while the corrections to reduction lead to a "background" of 0{1/N) which is 
approximately independent of loop size. As the coefficient of the decay — the potential, or 
more generally a hadron mass — increases, one has to go to ever higher values of A^. This 
problem should be less serious, however, for light particles, such as one expects to find as k 
approaches Kc- Indeed, calculating the pion mass would allow an important cross-check on 
our preferred interpretation that the system is confining, and chiral-symmetry breaking, in 
the infrared. 

One can also use the Nf = 2 AEK model away from the critical line as a single-site model 
whose long distance physics is that of the pure gauge theory. In other words, heavy adjoint 
fermions resolve the problems of the original Eguchi-Kawai model. The same holds true for 



the Nf = 1 theory [sl, 



To address the unresolved issues described above, one will need either to work at larger A^ 
or move to models with more than one site. The latter option seems most practical, and also 
has the advantage of being simpler to parallelize. First steps in this direction have been taken 
in Ref. It may also be advantageous to use twisted boundary conditions, as has been done 



for the Nf = 1 theory in Ref. since these appear to reduce the power of the corrections 



from 1/A^ to l/N"^. One can also consider using improved gauge and fermion actions, since 



these are known to clarify the infrared behavior in large-volume simulations 
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Appendix A: Models for eigenvalues of Dw 

In this appendix we describe various possible behaviors of the link eigenvalues and their 
implications for the eigenvalues of D]y. These models guide the interpretation of the results 
presented in Sec. IV Al for the spectrum of Dw. 

We consider the extreme weak coupling limit, 6 ^ 1, in which we must choose links that 
maximize S'gauge- This is achieved by links which can simultaneously diagonalized by a gauge 
transformation, i.e. for which one can have 

J9}. J9?. jeff' 



Uf, = diag[e'"^,e'"^,...e'"^ j V/i. (Al) 

What is needed for Dyy is the link in the adjoint representation [see Eq. (12. 4p ]. It is convenient 
to add a singlet and consider the link in the reducible N ^ N representation. In this case, it 
has composite indices, A = {ai, 02), with aj = 1, N, and is also diagonal: 

(^r"),, = (^.)... ^ f^r" = diag (. . . , e^^^-'^\ ...). (A2) 

Inserting this into the massless Wilson-Dirac operator, one findj 



4D^(mo = 0) = diag . . , | (4 - ^ cos 9;^'^') + ^ E ^T" | ' 



(A3) 



We multiply by 4 = 1/[2kc(3^ = 0)] in order that (in large volume) the operator becomes the Dirac operator 
with standard normalization in the naive continuum limit. This undoes the standard renormalization of the 
fermion fields by \/2k that is used to write D\y in the form of Eq. (|2.4p . 



46 



where we are using the abbreviation 

= e;^ - e;- . (A4) 

Thus D\Y is diagonal in color space, but not in Dirac space. The eigenvalues of D\Y{mQ = 0) 
are as follows. There are AN zero-modes, occurring when oi = as so that 6'^i"2 _ q_ Four 
of these are from the singlet, which we can now remove, leaving 4(A^ — 1) from the adjoint. 
The remaining 4N{N — 1) are each doubly degenerate (due to charge conjugation symmetry), 
come in complex conjugate pairs (due to 75-hermiticity), and have values 



Ka, = (4 - 5^ COS 6';^'^^) ±^ Yl ^i^' ^2'"' ' ^ "2) . (A5) 

The form of in Eq. f lA3p is exactly that of a four-dimensional free massless Wilson-Dirac 
fermion, with the momenta in lattice units ap^ replaced by 0^^"^. This is the standard way 
in which the large volume appears in the weak coupling limit, i.e. link eigenvalue differences 
become momenta If these eigenvalue differences are distributed such that the resulting 
"momenta" lie uniformly throughout a four-dimensional Brillouin zone, then the single-site 
Dw will approximate that of large-volume four- dimensional theory. The resulting spectrum 
has the five fingers mentioned above, whose "tips" occur when all four ap^ equal or tt. The 
tips are distinguished by the number of ap^ which equal — either 0, 1, 2, 3 or 4. Alternatively, 
if the eigenvalues are correlated in some way, then the spectrum will not, in general, have five 
fingers. There would, for example, be fewer fingers if the effective dimensionality is less than 
fourl^^ 

To illustrate these comments we discuss the results from three simple (and somewhat 
artificial) models for the eigenvalue distributions. In the first, we choose the eigenvalues in 
each direction to be evenly spaced around the unit circle, but in a randomly permuted order, 
with the permutation being independent in each direction. This leads to 

^r' = ^[^.(«i)-^MM], (A6) 



Given the weak-coupling form (|A5p . the spectrum is necessarily confined to lie between the ellipse {R — 
4)2 + (2/)2 = 16 and the four circles {R - i?o)^ + P = I, with Rq = 1, 3, 5, and 7. Here R and / are 
respectively the real and imaginary parts of the eigenvalues. Thus the distribution is kinematically forced 
to lie in one of the five fingers once |/| < 1. The presence of fingers per se is thus not significant, but the 
number of fingers which are populated is significant. 
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with cr^ a permutation of 1 — A^. In the eigenvalue-momentum correspondence these are the 
subset of the momenta available on an N"^ lattice with periodic boundary conditions. The 
Brillouin zone of such a lattice contains A^^ momenta, which is much larger than the A^^ 
values of 6*^^"^ produced by a single configuration. We thus assume further that independent 
configurations lead to independent permutations. Then, with of 0{N'^) configurations, one 
obtains the spectrum shown in Fig. [22] by the large (blue) dots. This is compared in the 
figure to the full spectrum of a 4-d free Wilson-Dirac operator on an A^^ lattice. One sees 
the appearance of the desired 5 fingers, but also that some points in the full spectrum are 
missing. This is due to the fact that ^^i"^ _ q only if ai = 02, in which case 0'^^°'^ = for 
all fi. The model thus cannot produce momenta proportional to (0,0,0,71,4), (0, 0, 71,3, 71,4) or 
(0, n2, 71,3, 71,4) (or their permutations). These "missing modes" have the largest impact on the 
left-most (and thus physical) finger, but become increasingly unimportant as A^ increases. We 
note that the distance of the fingertips to the real axis scales as 1/A^. 

In this model center symmetry is unbroken, and in particular all traces of open loops, Kn, 
vanish (unless all four are integer multiples of N). It gives an example where the effective size 
is Leff = A^, in the sense that the fermion operator after averaging over configurations has the 
same spectrum as a theory with volume A^^. The model is artificial in that eigenvalues do not 
fall on the "clock" values in unbroken phase, but, as seen above, are spread nearly uniformly. 
Nevertheless, it indicates how the spectrum of Dw can teach us about the distribution of 
eigenvalues and the effective dimensionality. 

Our second model is a variant of the first in which the eigenvalues still take clock values, 
but they are fully correlated between the directions. This breaks the Zj^ center symmetry 
down to the "diagonal" subgroup. More precisely, we assume that 6'^ = 6*" + {2tc / N)ni,^^ 
for all a, /i and z/, with n,^^^ an integer. This leads to 

= f (A7) 

where now there is a single permutation a for all four directions. This is the type of "locking" 
found in the quenched EK (QEK) model Q]. It leads to the spectrum of being that of 
a 1-d free Wilson fermion in a periodic box of length Lgg = A^, scaled up by a factor of 4. In 
this case the central three fingers are missing. 
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FIG. 22. Spectrum of 4L>vi/(0) in the model given by Eq. (jA6l) for = 16 and with 0{N^) con- 
figurations (large [blue] dots) compared to the spectrum of the free Wilson-Dirac operator on a 16^ 
lattice (small [red] dots). Only eigenvalues with positive imaginary part are shown; the spectra are 
symmetric under reflection in the Re A axis. 



The third model is inspired by the analysis of Ref. p^, in which it is shown that, at 
extremely weak coupling, the repulsion between eigenvalues leads to the formation of a four- 
dimensional crystal ii K = N^/'^ is an integer, and an approximately uniform distribution 
for other values of N . It is important to note that the analysis of Ref. 12|] holds only if 
the coupling h evaluated at the scale l/(aLeff) is much larger than unity. This requires that 
the lattice coupling, 6(l/a), grows logarithmically with N. This is not the standard limit in 
which reduction holds, in which h is fixed. Indeed, volume independence does not hold in this 
regime. We can express this realization of eigenvalues as 



\na na na nai 



mod(a — 1, i^),mod( ,K), mod( 



K 



^ i^),mod(^,i^) 



which breaks the center symmetry down to Zj^, but it is argued in Ref. 12| that fluctuations 
in the eigenvalues can lead to an averaging over the different "crystals" related by trans- 
formations, and thus the restoration of the full symmetry. The momenta 6'^i''^2 g^j^g ^j^^g same 
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for all crystals, and lead to a four- dimensional spectrum for -DvfUj The difference from the 
first model is that instead of the spectrum being that on a lattice of size L^s = it is on the 
much smaller size L^g = K = N^^^. This is the same scaling as for the orbifold construction 
discussed in the main text, and is the most conservative possibility. 

For the values of N that we use, one has L^s < 3. For such a small lattice, even though it 
is four-dimensional, the spectrum shows no fingers, i.e. no eigenvalues close to the real axis. 
Thus if this model of the eigenvalues provides even an approximate description of our data, 
we would not expect to see fingers. 

It is interesting that the above distribution of eigenvalues was first suggested in the context 



of the space embedding into color space of the Q 



For a further discussion on this point see Ref. 29| where two of us analyzed this eigenvalue 



IK model 



and then analyzed in Ref. 



distribution within the QEK model (where it is referred to as the 'Brillouin Zone' distribution). 
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